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ABSTRACT 



We have analyzed high-resolution spectra of relatively cool blue stragglers and main sequence 
turnoff stars in the old open cluster M67 (NGC 2682). We attempt to identify blue stragglers 
whose spectra are least contaminated by binary effects (contamination by a binary companion 
or absorption by circumstellar material). These "best" stragglers have metallicities ([Fe/H] = 
-0.05) and abundance ratios of the blue stragglers are not significantly different from those of the 
turnoff stars. Based on arguments from hydrodynamical models of stellar collisions, we assert 
that the current upper limits for the lithium abundances of all blue stragglers observed in M67 
(by us and others) are consistent with no mixing during the formation process, assuming pre- 
main sequence and main sequence depletion patterns observed for M67 main sequence stars. We 
discuss composition signatures that could more definitively distinguish between blue straggler 
formation mechanisms in open cluster stars. 

We confirm the spectroscopic detection of a binary companion to the straggler S 1082. From 
our spectra, we measure a projected rotational speed of 90 ± 20 km s _1 for the secondary, and 
find that its radial velocity varies with a peak-to-peak amplitude of > 25 km s — 1 . Because the 
radial velocities do not vary with a period corresponding to the partial eclipses in the system, we 
believe this system is currently undergoing mass transfer. In addition we present evidence that 
S 984 is a true blue straggler (and not an unresolved pair) . If this can be proven, our detection 
of lithium may indicate a collisional origin. 



Subject headings: stars: abundances — stars: blue stragglers — open clusters: individual (NGC2682, 
NGC7789) 



1. 



Introduction 



been identified in most, if not all, globular clusters 
that have accurate broad-band photometry. Blue 
straggler candidates have also been identified in a 
large number of open clusters, e.g. Ahumada & 
Lapasset (1995). 



Blue straggler (BS) stars were first identified 
as an unusual subclass of stars in the cluster M3 
(Sandage 1953). Since that time, stragglers have 



1 Based in-part on observations obtained with the 
Hobby-Eberly Telescope, which is a joint project of the Uni- 
versity of Texas at Austin, the Pennsylvania State Univer- 
sity, Stanford University, Ludwig-Maximillians-Universitat 
Miinchen, and Georg-August-Universitat Gottingen. 



Blue stragglers are commonly defined as those 
stars that are brighter and bluer (hotter) than 
the main sequence (MS) turnoff of the majority 
of the cluster stars, but along an apparent exten- 
sion of the main sequence. At present, the lead- 



ing explanation involves mass transfer in and/or 
merger of a binary star system, or the collision of 
stars (whether or not a binary is involved). These 
mechanisms can create a main sequence star with 
a mass greater than would be expected given the 
age of the cluster. 

The above definition of blue straggler is overly 
restrictive if this is the correct explanation. Un- 
doubtedly the same processes act on lower mass 
stars, although the remnant stars would not be 
easily identifiable if they were not brighter than 
the turnoff luminosity. The color restriction is 
also unrealistic since real scatter in the colors of 
blue stragglers in globular clusters has been ob- 
served, e.g. Fusi Pecci et al. (1992). Sills & Bai- 
lyn (1999) find that this color scatter results from 
a lack of large-scale mixing during the creation 
process, which means that the amount of helium 
in the core can vary from straggler to straggler. 

Hydrodynamical studies (Lombardi et al. 1995; 
Sandquist et al. 1997) earlier predicted that large- 
scale mixing does not occur during stellar colli- 
sions. A collisionally-produced straggler is not ex- 
pected to burn significant lithium initially even if 
high enough temperatures are reached because the 
timescale for the straggler to thermally adjust to 
a new equilibrium configuration is much shorter 
than the lithium destruction timescale Sills et al. 
(1997). Binary mass transfer is likely to lead to 
a more easily observable abundance pattern since 
the last of the transferred gas should have come 
from deep within the parent star (e.g, Sarna & de 
Greve 1996). 

In this paper we have attempted to examine 
this question by taking spectra of blue straggler 
stars in the open cluster M67 (NGC 2682). M67 
makes an ideal target because it is relatively close 
(V-M v w 9.7; Dinescu et al. 1995), and also fairly 
old (~ 4.0 Gyr; Dinescu ct al. 1995). In addition, 
the density of stars in open clusters is lower than 
in nearly all globular clusters, making binary evo- 
lution and collisions involving binaries the most 
likely formation mechanisms (Leonard & Linnell 
1992). Globular cluster blue stragglers are not yet 
practically observable because of their distances. 
Although much attention has been spent on M67 
blue stragglers, few observations are available for 
the coolest ones. 

Few blue stragglers have abundance analyses. 
Mathys (1991) observed a sample of 11 blue strag- 



glers in M67, although only 2 were examined in 
detail. Andricvsky (1998) observed 4 blue strag- 
glers in Praesepe (NGC 2632) from the catalog of 
Ahumada & Lapasset (1995), although only one 
(HD 73666) is undeniably a straggler according to 
its photometry. In observing the blue stragglers 
of M67 we are most interested in species such as 
Li, C, N, and O that are processed by nuclear re- 
actions at different temperatures (and therefore, 
depths) in a star. The reddest blue stragglers in 
the cluster have effective temperatures that can 
allow us to get abundances for all of these ele- 
ments, and are close enough to the temperatures of 
turnoff stars that we can straightforwardly make 
comparisons with stars having (presumably) un- 
mixed envelopes. 

In §2 we describe the spectroscopic observa- 
tions, and in §3 we present the abundance analysis 
for our sample. In §4 we discuss the constraints we 
can place on blue straggler formation mechanisms 
in this cluster. 

2. Observations 

Table 1 displays the program stars, basic pho- 
tometry, and the signal-to-noise obtained. The 
identification numbers we use are from the proper 
motion study of Sanders (1977). We chose our 
blue straggler candidates from Ahumada & Lapas- 
set (1995) and our MS members from the proper 
motion study of Zhao et al. (1993). All of the cho- 
sen M 67 stars are 99% probability proper motion 
members. The photometry for our program stars 
was taken from Fan et al. (1996), Montgomery et 
al. (1993), Chevalier & Ilovaisky (1991), Gilliland 
ct al. (1991), and Sanders (1989). The B — V col- 
ors of Sanders and Gilliland et al. were adjusted 
by 0.006 and —0.002, respectively. The V mag- 
nitudes of Sanders, Gilliland et al. and Cheva- 
lier & Ilovaisky were adjusted by 0.017, —0.009, 
0.012, respectively. Several measurements made 
by Sanders were found to be in disagreement with 
the other sources, including ancillary measure- 
ments by Anupama et al. (1994), Murray et al. 
(1965), and Murray & Clements (1968). For this 
reason, the Sanders measurements for S 1183, S 
821, S 975, and S 984 were not used to derive 
the average photometry. The final adopted values 
are given in Table 1, and the CMD of the cluster 
and the observed blue stragglers are shown in Fig- 
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ure 1. We observed one blue straggler candidate 
(Ahumada & Lapasset 1995) in NGC 7789. Due 
to poor weather additional NGC 7789 candidates 
were not observed. 

The majority of the observations were taken 
with the McDonald Observatory 2.7 m "2d-coude" 
(Tull ct al. 1995). The resultant spectra have 
R = 30, 000 and cover from 3800 A to 10100 A. 
Each star was observed on several different nights 
in order to look for radial velocity variations. Ra- 
dial velocities were obtained from individual spec- 
tra, but the S/N quoted in Table 1 are those of 
the combined spectra. The combined spectra were 
also used for the abundance analysis. Additional 
observations were obtained with the Hobby-Eberly 
Telescope (Ramsey et al. 1994, HET) during com- 
missioning of the UFOE (Harlow et al. 1996) spec- 
trograph in early 1999 and in the first science pe- 
riod of 2000 with the UFOE on the HET. The 
HET spectra were used to obtain radial velocities 
and have lower resolution (R = 11,000) than the 
2.7 m spectra, and so are not included in the com- 
bined spectra. The HET spectra cover from 4600 
A to 9100 A. In March 2000, the McDonald Obser- 
vatory 2.1 m Cassegrain echelle was used to obtain 
a few additional observations of two of the pro- 
gram stars. Despite having resolution higher than 
the 2.7 m "2d-coude" spectra, these were not in- 
cluded in the abundance analysis because the slit 
length did not allow proper sky subtraction, which 
would not interfere with the radial velocity study 
since M67 is well separated from the solar rest ve- 
locity. The 2.1 m spectra cover from 5400 A to 
6600 A. Table 2 lists the heliocentric Julian dates 
of our observations. In addition to the program 
stars one additional M67 star (MMJ5554) was ob- 
served. This star is only 3" away from S 975 and 
its color makes it either a turnoff star or a pos- 
sible blue straggler candidate. A set of bright, 
rapidly-rotating hot stars were observed at vari- 
ous air masses with both telescopes for sky line 
identification. Figure 2 exhibits a sample of the 
spectra obtained on the McDonald 2.1m. 

3. Analysis 

3.1. Radial and Rotational Velocities 

The spectra were analyzed using the IRAF 
echelle package. Initial data reduction techniques 
included bias removal, flat field correction, scat- 



tered light removal, sky subtraction and wave- 
length calibration. To remove instrumental wave- 
length zero-point errors the spectra were shifted 
in velocity so that the telluric features fell at their 
rest velocities. The largest zero-point correction 
was 1.1 km s _1 . Each program star spectrum was 
divided by the hot star spectrum with the most 
similar telluric line depths. The spectra were then 
cross-correlated against a synthetic spectrum to 
obtain the radial velocities, which were then cor- 
rected to the heliocentric rest frame. A single syn- 
thetic spectrum was used for all the stars (T e ff — 
6400 K, log g = 4.3, v t = 2.0, [Fe/H] = -0.05, and 
no instrumental or rotational broadening). The 
code used to create the synthetic spectrum is the 
latest version of MOOG (Sneden 1973). To deter- 
mine the radial velocities of the broad component 
of S 1082 we took an additional step to remove 
the narrow component. We fit the spectrum with a 
very high order spline (high enough to remove any 
broad component) and then divided this normal- 
ized spectrum into the original spectrum leaving 
only the blaze function and any broad line compo- 
nents. These cleaned broad-line spectra were then 
subjected to the analysis described above. Table 2 
lists the heliocentric radial velocities and errors for 
each spectrum. The mean radial velocity of M67 is 
33.6 ± 0.72 km s" 1 (Girard et al. 1989). We detect 
3 binary periods not previously known for stars in 
our sample. S 821 and S 984 are found to have 
radial velocity variations with periods of 26.259 ± 
0.002d and 1.465 ± 0.001, respectively. Figure 3 
shows the phased radial velocity data from Math- 
ieu et al. (1986) and this work for these periods. 
S 821 probably has an eccentricity of 0.4 ± 0.1, 
while the data from S 984 is consistent with a cir- 
cular orbit. S 1082 is found to be a double lined 
spectroscopic binary and will be discussed in §3. 

The 2.7 m "2d-coude" spectra were combined 
after shifting to the heliocentric rest velocity. New 
synthetic spectra were created with T e ff = 6400 
K, log g = 4.3, v t = 2.0, [Fe/H] = -0.05, a Gaus- 
sian R = 30, 000 profile, and various rotational 
velocities. The unbroadened spectrum was cross- 
correlated against the broadened spectra and the 
full width half maxima (FWHM) of the correla- 
tion peaks were measured. These FWHM were 
then compared to the FWHM of the correlation 
peaks created by the cross-correlation of the com- 
bined program star spectra and the unbroadened 
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synthetic spectrum. From this comparison we de- 
termine the rotational velocity of the combined 
spectrum. We adopt a conservative lower limit (9 
km s _1 ) to our ability to measure rotational veloc- 
ities. Table 1 lists the rotational velocities found 
with this method. 

We detect rotation in our spectra for only 3 of 
the blue straggler stars in our sample: S 975, S 
997, and S 1082. All have relatively low rotation 
speeds. In addition, we marginally detect rotation 
in one of the main sequence stars (S 1271). The 
signal-to-noise ratios for the main sequence stars 
was comparable to those for the blue stragglers. 

We observed one blue straggler candidate in 
NGC 7789, M 1251, but found it to be a cluster 
non-member based on its radial velocity. 

3.2. Abundance Analysis 

The spectral lines used in our analysis were 
taken from Shctronc (1996), Mathys (1991), and 
Edvardsson et al. (1993). The oscillator strengths 
were taken from these sources as well as the 
National Institute of Standards and Technology 
Atomic Spectra Database. Each line's equivalent 
width was measured with the IRAF task splot. 
These lines, oscillator strengths, and equivalent 
widths are listed in Table 3. We employed the 
latest version of the LTE spectral analysis code 
MOOG (Sneden 1973) and the Kurucz (1993) grid 
of ATLAS models in an iterative abundance anal- 
ysis. 

The initial estimate of the effective temperature 
for each star was determined from the dereddened 
B — V color (we assumed E(B - V) = 0.05). Ef- 
fective temperatures were then determined from 
the relationship T eff = 1808(5 - Vf - 6103(B - 
V) + 8899 (Sodcrblom et al. 1993). Each effec- 
tive temperature was then fine tuned by forcing 
the slope of abundances from Fc I lines versus ex- 
citation potential to be zero. This tuning of the 
temperature was restricted to the range of colors 
consistent with the errors in the photometry. The 
two largest corrections were found for S 1082 (152 
K; likely binary blending of colors, see §3) and S 
815 (42 K). The surface gravity was set by en- 
forcing abundance equilibrium between the Fe I 
and Fe II lines. Based upon the errors in the Fe 
I and Fe II abundances we estimate the average 
uncertainty in the surface gravity to be 0.15 dex 



(1 sigma). The microturbulent velocity was deter- 
mined by forcing the slope of Fe I line abundances 
versus equivalent width to be zero. These adjust- 
ments to the spectral model for each star were 
done iteratively until all the parameters (T e ff, log 
g, v t ) fulfilled the above requirements within the 
abundance uncertainty of each spectrum. The de- 
rived atmospheric parameters and abundances for 
each star are listed in Table 4. 

Synthetic spectra were also employed to verify 
the abundances of the C I lines. Using the derived 
atmospheric parameters for each star a synthetic 
spectrum was generated and plotted over the ob- 
served spectrum and visually inspected. This 
method did not yield any significant abundance 
differences. 

The three O I lines found at 7773A are known 
to need NLTE corrections. These corrections can 
become very significant at high temperatures (T 
> 6800 K). Over the narrow range of tempera- 
tures we considered the NLTE corrections should 
be relatively small. However, we have made some 
attempt to correct for NLTE effects using the 
lookup table employed by Gratton et al. (1999). 
Even if the magnitude of our corrections are incor- 
rect, the relative oxygen abundances of the stars 
can be viewed with some confidence. In Figure 4 
we display the combined equivalent widths of the 
O I lines plotted against their colors. The stars 
from Mathys (1991) are shown as squares. The 
large crosses are taken from the Varenne & Monicr 
(1999) study of oxygen in the Hyades. Varenne & 
Monier (1999) found a constant abundance along 
the main sequence. Note the rapid rise in O I EW 
near B — V = 0.43 or T eff = 6840 K, caused by 
the onset of NLTE effects. The open and filled cir- 
cles represent abundance analyses with constant 
oxygen abundance ([O/Fe] = 0.0) for several tem- 
peratures with the NLTE corrections of Gratton 
et al. (1999) and Faraggiana et al. (1988), respec- 
tively. A summary of our oxygen abundances is 
given in Table 4. 

4. Discussion 

Because the stragglers in M67 have been ob- 
served extensively, we have a better idea of their 
present status than the stragglers in any other 
cluster. It is therefore worth briefly summarizing 
what is known about the stars in our sample so 
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that we can test whether the compositions reflect 
differences in straggler formation mechanisms. 

4.1. Discussion of Individual Stars 

4-1.1. Blue Stragglers 

S 975 (F 90): This star is known to be a 
binary (P = 1221 d) with small eccentricity 
(e = 0.088 ± 0.060) from the work of Latham 
& Milone (1996). (Leonard 1996) indicates that 
the system is likely to have resulted from case C 
(asymptotic giant branch) mass transfer. Lands- 
man et al. (1998) assert that there is strong evi- 
dence of a hot companion in their ultraviolet pho- 
tometry, and conclude that the system is likely to 
be the result of Algol-type binary mass transfer, 
although they adopt a photometric temperature 
(Mathys 1991) that is more than 300 K cooler 
than our derived value. Mathys states that the 
photometry could be contaminated by a nearby 
star, so that his temperature estimate is question- 
able. The contaminating star must be cooler than 
S 975 in order to produce the photometric tem- 
perature. MMJ 5554 (V = 12.839, B — V = 0.532) 
is probably the cause since it is only 3" away from 
S 975, and its colors place it very near the bluest 
extension of the turnoff. Because our spectral tem- 
perature did not need to be significantly corrected 
from our BV photometric temperature, it is un- 
likely that this star contaminated our spectrum. 
Judging from Figure 2 of Landsman ct al. (1998), 
it probably did not contribute significantly to the 
UIT UV flux. 

According to Landsman et al., "uncertainty of 
200 K in the value of T e f f corresponds to about a 
0.5 mag uncertainty in the predicted UIT magni- 
tude" , and "the UIT magnitude of F90 is proba- 
bly uncertain by close to a factor of 2" because it 
was contaminated by the light of the bright blue 
straggler S 977. Based on this and our higher 
T e ff, the evidence for a hot companion to S 975 
is considerably weaker, although there is still an 
unaccounted-for UV excess. The agreement be- 
tween our photometric and spectroscopic temper- 
atures indicates that the companion does not con- 
tribute significantly to the visual flux (consistent 
with being a white dwarf). However, we have been 
unable to detect the companion's spectral signa- 
ture, so we are unable to judge whether the spec- 
trum is contaminated. The star does have an un- 



usually large O I EW though (see §4.2.2), which 
may betray the presence of circumstellar material. 

S 984 (F 134): Because S 984 lies within 0.75 
mag of the fiducial line of the cluster in the color- 
magnitude diagram, it is possible that it is a non- 
interacting binary composed of two main sequence 
stars. The velocity dispersion in the radial veloc- 
ity data of Mathieu et al. (1986) was marginally 
higher than average. After combining the velocity 
data from Mathieu et al. (1986) with our derived 
velocities, we conducted a periodicity analysis and 
found a likely period of 1.465 ± 0.001 days. The 
peak-to-peak variation is between about 3 and 5 
km s _1 , depending on whether we accept one out- 
lying observation from Mathieu et al. 

The period, radial velocity variation, and pho- 
tometry do constrain the identity of S 984 though. 
If we hypothesize that S 984 is a binary composed 
of stars that have not interacted, it is possible to 
match the photometry with two main sequence 
stars having V « 12.7 and 13.4. This is the max- 
imum brightness contrast that is possible if the 
primary is to be on the fiducial line of the clus- 
ter. However, this constrains the mass ratio to be 
0.9 < q < 1.0. The orbital inclination is then con- 
strained to be within 1.5 degrees of face-on. The 
lithium abundance, which is consistent with that 
of a single turnoff star, is low for a tidally-locked 
binary (as is implied if the radial velocity varia- 
tions reflect orbital motions). So, there is a fairly 
strong (although circumstantial) case that S 984 is 
a true blue straggler, and not a photometric blend 
of two stars. There may still be a companion, but 
it is unlikely to contribute to the flux. If this is the 
case, S 984 is probably the result of a collisional 
merger, since binary mass transfer should result 
in a lack of surface lithium. For these reasons, we 
regard it as likely that the spectrum of S 984 is at 
most minimally contaminated by another star. 

S 997 (F 124): This star is known to be in an 
eccentric (e = 0.342 ±0.082) binary of period 4913 
d (Latham & Milone 1996). The eccentricity of 
the system probably rules out a history as a mass 
transfer binary. The fact that the primary is a 
blue straggler casts doubt on a binary merger since 
it would either have required the system to have 
been in a rare triple system, or to have tidally cap- 
tured a star after merger (Leonard 1996). Either 
way, we do not have evidence to judge whether its 
spectrum is contaminated. 
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S 1082 (F 131; ES Cancrii): Pritchet & 
Glaspey (1991) and Mathys (1991) reported that 
the spectrum of S 1082 has a possible composite 
nature. X-ray emission has been detected from 
this star by Belloni et al. (1993). Landsman et al. 

(1998) find evidence of a hot subluminous com- 
panion in their ultraviolet photometry. However, 
Landsman et al. (1998) base their predicted UIT 
magnitudes on a temperature 120 K cooler than 
found in our analysis. Since a 200 K change in 
temperature corresponds to 0.5 mi52 magnitudes, 
their evidence for a hot subluminous companion 
is weakened, although not eliminated. The dis- 
covery of a variable secondary Ha spectral fea- 
ture by Mathys (1991) and van den Berg et al. 

(1999) suggests the presence of a hot rapidly ro- 
tating secondary or an outflowing wind. Further, 
the spectroscopy of van den Berg et al. (1999) ap- 
pears to rule out magnetic activity as the source 
of the X-ray emission. This evidence supports the 
hypothesis that this is an Algol-type binary cur- 
rently undergoing mass transfer. 

In addition to being able to detect the second 
spectral component in Ha, the Na D lines, and the 
O I triplet as Mathys (1991) did, we also detected 
the secondary in H/3, H7, Paschen(8862), the Ca 
II IR triplet, lines of moderately strong neutral el- 
ement species (e.g. Fe I, Ca I) and a few strong Fe 
II lines. The abundances derived for 1082 are 0.2 
dex lower than the cluster mean. This can be ex- 
plained if the second component produces 40-50% 
of the continuum flux but only a tiny (<~10%) con- 
tribution to the line flux. The rotational velocity 
implied for the second spectral feature is 90±20 
km s _1 , which explains the small contribution to 
the line flux. 

For the first time, we are able to see radial 
velocity variations in the secondary (broad-lined) 
component amounting to at least 25 km s" 1 . Our 
phase coverage is not good enough to state more 
than a lower limit. Figure 2 exhibits the scries 
of spectra we obtained of S 1082 on the McDon- 
ald 2.1m. The broad component can clearly be 
seen when compared with the spectra of S 984. In 
addition, the broad component of S 1082 moves 
redward over the 5 hours the star was observed. 
For the primary, we find a peak-to-peak amplitude 
of about 5.5 km s _1 . From the data of (Mathieu 
et al. 1986), there is a peak-to-pcak amplitude of 
about 7 km s^ 1 . 



Goranskij et al. (1992) found from photometry 
that the system has two partial eclipses per orbit 
of different depth, with a period of 1.0677978 d. 
Using this period (or half of it), we see no sign of a 
repeating pattern in the radial velocity data. Al- 
lowing the period to vary, we find a highest prob- 
ability period of 1.87 ± 0.05 d, with less likely 
values of 0.7 d and 2.3 d (see Figure 3). How- 
ever, none of these periods unveils a satisfactory 
pattern in the radial velocity data. We have two 
sequences of spectra taken on succeeding nights 
(centered around phase 0.5 of the photometric pe- 
riod) that indicate sudden (« 2 h) jumps in radial 
velocity for both components. Thus, the radial ve- 
locity data imply that we are seeing the signature 
of material outside of the photospheres of the two 
stars. 

The spectral temperature derived for the nar- 
row component is more than 100 K hotter than the 
photometry predicts. Using the color as the com- 
posite temperature for the system, our spectro- 
scopic temperature suggests that the narrow line 
component may be 200 K hotter than the broad 
lined component. Further evidence for a tem- 
perature difference between the two components 
comes from the spectrum: the broad component's 
high excitation lines of O I have smaller equiva- 
lent widths than those of the narrow component, 
while the low excitation Na D lines have higher 
equivalent widths. This implies a cooler tempera- 
ture, although strong Fe II lines may indicate that 
a lower surface gravity also comes into play. The 
relatively large brightness (My ~ 2.5) and hot 
temperature (T ~ 6800 K) of the secondary indi- 
cates that it has probably been severely disturbed 
by its interaction with the primary. 

S 2204 (F 130): We observed this star initially 
believing it to be a main sequence star. Evidence 
indicates that it is very probably a faint blue strag- 
gler, and it is identified as such in the catalog of 
Ahumada & Lapasset (1995). The proper motion 
studies of Girard et al. (1989) and Sanders (1977) 
cite membership probabilities of 98% and 83% re- 
spectively, and radial velocity measurements by 
Mathieu et al. (1986) and us put it right on the 
cluster average. The photometry of Montgomery 
ct al. (1993) puts it (B - V = 0.448) over 0.07 
mag to the blue of the bluest stars at the turnoff, 
and over 0.1 mag from the fiducial line of the clus- 
ter. The quoted mean error in the photometry was 
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0.016. Significant contamination of the spectrum 
by a main sequence companion is not likely since 
this would either make the star brighter than the 
turnoff, or closer in color. 

A hot, faint companion could also potentially 
create the observed colors. Although S 2204 was 
in the field of the ultraviolet observations of Lands- 
man et al. (1998), it was not detected. Their de- 
tection of several white dwarf candidates thereby 
places an upper limit on the temperature of a 
possible white dwarf companion. Landsman ct 
al. (1998) estimate their detection threshold for 
a white dwarf with a pure hydrogen atmosphere 
would be about 21,300 K. A white dwarf below 
the detection threshold would change the broad- 
band color by a negligible amount. 

For stars close to the turnoff, the possibility of 
delayed star formation (a very unlikely possibility 
for brighter blue stragglers) should be considered. 
For S 2204, a delay of more than 1 Gyr (about 25% 
of the cluster's age) would be necessary according 
to the isochrone fits of Dinescu et al. (1995) if S 
2204 is an undisturbed main sequence star. On 
this basis we discard the delayed star formation 
picture for S 2204. 

With the elimination of these hypotheses, the 
most likely explanation of S 2204 is as a blue strag- 
gler of the merged-star or mass-transfer variety. 
Although there is currently not enough radial ve- 
locity information to more strongly test for the 
presence of a companion to S 2204, we argue that 
the primary would be identified as a blue strag- 
gler even without one, and that such a companion 
would not significantly modify the spectrum of the 
straggler. 

4.1.2. Turnoff Stars 

S 821: We have detected binary motions for 
this object with a period P = 26.259 ±0.002 d and 
eccentricity e « 0.4±0.1 (Figure 3). Broad-band 
photometry places it very close to the fiducial line 
of the cluster. This probably indicates that the 
spectrum is uncontaminated by the companion, 
although there is a chance that it is composed of 
two fainter, nearly equal mass stars. 

4.2. Composition Analysis 

The blue stragglers we have observed can be di- 
vided into two probable categories: mass transfer 



binaries (S 975, S 1082), and collision products (S 
984, S 997) . At present we do not have enough in- 
formation to classify S 2204. We also have spectra 
for 5 stars near the main sequence turnoff (S 815, 
S 821, S 1183, and S 1271). Our goal is to com- 
pare the abundances to see if we can find chemical 
signatures of the blue straggler formation process. 

4-2.1. Lithium 

If significant mixing of a star's envelope has 
occurred, then the lithium abundance may show 
the results first since it is consumed at relatively 
low temperatures (T > 2.5 x 10 6 K). Jones et 
al. (1999) present the most recent compilation of 
lithium abundances for M67 stars. 

Our results are given in Table 4. Two of our 
stars (S 997 and S 2204) were previously observed 
by Garcia Lopez et al. (1988), who found upper 
limits. We derive upper limits that are 0.5 and 
0.8 dex lower, respectively. With the exception 
of S 984, we are only able to derive upper lim- 
its for the stars in our sample. The values for 
the main sequence stars are consistent with the 
lower envelope of lithium abundances (see Jones 
et al. 1999, for references to other papers on M67 
lithium abundances). 

Hobbs & Mathicu (1991) previously observed 
the relatively cool blue stragglers S 1072 and S 
1082, while Pritchet & Glaspey (1991) observed 
stragglers S 752, S 997, S 1082, S 1263, S 1267, S 
1280, and S 1284. Both sets of authors found only 
upper limits for their samples. So, the lithium 
abundances do not as yet provide a dependable 
way of distinguishing among different formation 
mechanisms. 

4.2.2. CNO elements 

The CNO elements can potentially constrain 
the amount of mixing that occurs in a star, since 
the CNO cycle starts to shuffle abundances at 
around 10 7 K. Surface material is unlikely to reach 
this temperature in main sequence stars in the 
mass range present in M67, even if violent inter- 
actions occurred. However, if mass transfer were 
to occur on the giant branch, processed material 
might be dropped onto the surface of the straggler. 

Our results for carbon and nitrogen are given 
in Table 4, while our oxygen results are found in 
Table 5. Our nitrogen abundances are very uncer- 
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tain because of the very small equivalent widths 
found at this temperature, so we will not dis- 
cuss them further except to note that the nitro- 
gen abundances for the blue stragglers appear to 
be lower (0.2 - 0.3 dex) than the abundances found 
in the turnoff stars. The carbon abundances for 
all of the stars in our sample fall very close to the 
solar ratio. Comparing the oxygen abundances 
of the blue stragglers as a group to those of the 
main sequence stars, the blue stragglers (dropping 
S 975) also have an average abundance near that 
of the turnoff stars. However, several stragglers 
have anomalously high or low values. 

The straggler S 975 is problematic to analyze. 
On one hand, the evidence points to mass trans- 
fer of the type that could produce an observable 
chemical signature. On the other hand, the com- 
posite nature of the system appears to be affect- 
ing our equivalent width measurements. The mea- 
sured O abundance is rather high compared to all 
of the other stars in our sample. The expectation 
is that the mass transfer should lead to a deple- 
tion of oxygen if there was CNO processing. The 
anomalous abundance is possibly due to excess line 
flux contributed by a companion star or to inaccu- 
rate NLTE abundance corrections since S 975 lies 
very near the point where the NLTE corrections 
become large. 

For the two stars (S 975 and S 1082) that 
overlap with the blue straggler sample of Mathys 
(1991), the agreement of the O I triplet equiva- 
lent widths is very good. We have attempted to 
test the hypothesis that the blue stragglers and 
main sequence stars observed all have the same 
oxygen abundance. As seen in Figure 4, the ma- 
jority of the stars have EWs that are consistent. 
The exceptions are S 975 (which appears to have a 
companion that affects the line flux), S 1434 (B-V 
= 0.12, EW(0 I) = 1071 Mathys (1991)) (which 
also has a high EW, possibly related to a post- 
main sequence companion inferred from infrared 
photometry; Peterson et al. 1984), and S 968 (B- 
V = 0.12, EW(0 I) = 687 Mathys (1991)). S 968 
is known to be an Am star, which may account for 
the difference in that case. However, S 752 (B-V 
= 0.29, EW(0 I) = 754 Mathys (1991)) is also 
an Am star, and it does not appear to have an 
anomalous equivalent width. Finally, we find that 
S 1082 has a low O EW because of the continuum 
flux contribution of an nearly equal brightness sec- 



ondary. 

To aid in the interpretation of this figure we 
compare stars from other studies. The large 
crosses are taken from the Varenne & Monicr 
(1999) study of oxygen along the main sequence 
of the Hyades. Varenne & Monier found a nearly 
constant abundance along the main sequence. 
Note the rapid rise in O I EW at B - V = 0.43 
{{B - V)o = 0.38) or T eff = 6840 K due to 
NLTE effects. The open circles and filled circles 
represent synthetic abundance analysis employ- 
ing constant oxygen abundance ([O/Fe] = 0.0) for 
several temperature with the NLTE corrections of 
Gratton et al. (1999) and Faraggiana et al. (1988), 
respectively. 

To the accuracy of our measurements, we see no 
evidence for differences between the abundances of 
the blue stragglers and the turnoff stars. 

In an analysis of two hotter blue stragglers in 
M67, Mathys (1991) found total CNO abundances 
that were lower than those of giants in the cluster. 
Carbon and oxygen in particular were significantly 
depleted relative to iron. None of the blue strag- 
glers in our sample show this pattern of depletion. 
Nor do our blue stragglers show differences in C 
and O abundances compared to our turnoff stars. 
The abundance results are thus consistent with a 
collisional formation. However, the low eccentric- 
ities of some of the binaries containing stragglers 
are difficult to explain in that case. 

4.2.3. Other Elements 

Nuclear processing near the hydrogen burning 
shell of a giant star could produce anomalous 
abundances in other, more easily observable ele- 
ments. Surface anomalies in globular cluster giant 
stars have been observed in a number of globular 
clusters (e.g. M3, M13, M15, M71, M92, 47 Tuc). 

Gilroy & Brown (1991) showed that the carbon 
isotope ratio C 12 /C 13 in M67 giants is lower than 
can be predicted by standard evolutionary theory. 
However, the required extra mixing does not ex- 
tend into regions where the oxygen abundance is 
modified, as is found in metal-poor globular clus- 
ters. In fact, the amount of mixing seen in disk 
giants (Lambert & Ries 1981) is small in compari- 
son to globular cluster giants (Shetrone 1996). Na, 
an element which is burned at a relatively low tem- 
perature (Langer et al. 1997; Cavallo et al. 1998), 
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has not been studied well in open cluster giants. 

Here in particular our ability to discern which 
stragglers arc most likely to be free from con- 
tamination is important. The best candidates (S 
984 and S 2204) show fairly minor differences in 
abundance with the turnoff stars. We must also 
be aware that there may also be selection effects 
in that stragglers that form by some mechanisms 
may have unavoidably contaminated spectra. The 
straggler S 1082 shows abundances that set it off 
from the other stars in our sample — in particular 
Ni and Ba have lower than expected abundances. 
The measured iron abundance is very different 
from all of the other cluster stars. The tempta- 
tion is to think that this may relate to the mass 
transfer that is inferred. But as mentioned above 
this star is part of a binary having a composite 
spectrum that is difficult to disentangle. 

4.3. Interpretation 

4- 3.1. Lithium 

To date there has not been a detection of 
lithium in a blue straggler in an open cluster. 
There are several possible reasons for this, some 
of which have interesting consequences. 

If two stars were to mix completely during a 
collision, depletions of a factor of 50 might be ex- 
pected (Pritchet & Glaspey 1991) because lithium 
would get diluted throughout the blue straggler. 
However, hydrodynamical studies of stellar colli- 
sions (Lombardi et al. 1995; Sandquist et al. 1997) 
indicate that the collision itself is unlikely to mix 
the stars substantially — the high entropy gas near 
the surfaces of the input stars remains at high 
entropy near the surface of the remnant (Lom- 
bardi et al. 1995). As a result, a collisionally pro- 
duced straggler is not expected to burn significant 
lithium initially, even if high enough temperatures 
are reached (Sills et al. 1997). 

For blue stragglers at or above the turnoff of 
the cluster that were formed via stellar collisions, 
the lithium abundance observed at the surface is 
therefore likely to be the same as the abundance 
at the surface of the input star with the higher 
surface entropy. For main sequence stars, this 
would be the more massive star. If the difference 
in mass of the input stars is relatively large, gas 
from the lower mass star is likely to be completely 
sequestered in the interior of the straggler where 



it would be unobservablc. 

For open clusters, it is well known that lithium 
abundance drops with decreasing effective temper- 
ature along the main sequence (e.g. Jones et al. 
1999). For an open cluster as old as M67, two- 
body collisions that produce an identifiable blue 
straggler could involve a star of as low a mass as 
0.6 M as the more massive star. According to 
standard stellar models, main sequence stars with 
masses less than about 0.8 M Q are likely to have 
burned enough of their lithium that it would be 
undetectable. Standard models are currently un- 
able reproduce the lithium abundance patterns in 
the Hyades though. Extrapolating from Hyades 
measurements, stars with masses less than about 
1.0 M Q are likely to be undetectable. 

While collisions of 0.3 M and a 1.0 M Q 
stars could create blue stragglers with observable 
lithium, two factors work against this. First, such 
collisions are relatively rare thanks to the bias in 
the cluster luminosity function toward stars rela- 
tively near the turnoff (Montgomery et al. 1993). 
Second, significant main sequence Li depletion is 
possible after the blue straggler is formed. In the 
Hyades, a dip in lithium abundances is observed 
for stars with T e fj w 6600 K (the total range is 
about 300 K). It is believed that this feature may 
be due to meridional circulation or shear turbu- 
lence (for a recent review, see Talon & Charbonnel 
1998). For M67 no current main sequence mem- 
bers fall in this region, but some blue stragglers 
do. S 997 and S 2204 fall in the Li gap and 
may have depleted their Li even if they arrived 
on the blue straggler sequence with a significant 
supply. The age of the Hyades (~ 625 Myr; Per- 
ryman et al. 1998) provides an upper limit for 
the lithium destruction timescale for these stars if 
this mechanism does act. The lifetime of a blue 
straggler falling in that range of effective temper- 
atures would be about 3 Gyr, in which case there 
is less than 25% probability of finding lithium. 
Blue stragglers which form with more helium in 
their cores (because the lower mass input star was 
relatively evolved) would have higher probabilities 
of showing lithium, but shorter lifetimes as blue 
stragglers. 

More massive blue stragglers (T e ff > 7000 K) 
would be more likely to retain surface lithium be- 
cause both the straggler and the input stars lack 
extensive surface convection zones, and because 
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the stragglers are hotter than the lithium gap. A 
possible problem is the mass loss associated with 
a collision. In hydrodynamical models of globular 
cluster star collisions by Lombardi et al. (1995) 
and Sandquist et al. (1997), this mass loss was less 
than approximately 6% of the total input mass, 
with the largest mass loss occuring for relatively 
rare head-on collisions. However, mass loss would 
tend to remove the most lithium-rich material (for 
example, see Table 4 of Lombardi et al. (1995)). 
For stars with masses between about 1.0 and 1.25 
M Q (approximately the turnoff mass), the lithium- 
rich regions of the star add up to approximately 
4-5% of the mass. 

Blue stragglers formed via binary mass trans- 
fer are unlikely to show surface lithium (Sills et 
al. 1997). In principle, the brighter, more massive 
blue stragglers are the ones to observe to look for 
lithium evidence that different mechanisms are re- 
sponsible for the formation of M67 blue stragglers. 
Unfortunately, the Li line is very temperature sen- 
sitive and becomes extremely weak in hot massive 
stars. 

4.3.2. CNO elements 

Although they are less likely to show the effects 
of relatively low levels of mixing, CNO elements 
may prove to be more useful in distinguishing be- 
tween stragglers that have formed by collisions and 
by binary mass transfer. Because the CNO ele- 
ments undergo nuclear processing at higher tem- 
peratures than lithium, it is very unlikely that the 
surface abundances will be modified as the result 
of a collisional merger. Massive input stars will 
not have modified their surface CNO abundances 
because their convection zones are not extensive 
enough, while low mass input stars will in general 
become completely engulfed in a higher mass star 
due to their lower entropies. 

On the other hand, if Algol-type binary mass 
transfer has occurred (on the red giant branch or 
asymptotic giant branch), the mass visible at the 
surface of the blue straggler is taken from deep 
inside the primary where it could have undergone 
nuclear processing. Spectroscopy of giant branch 
stars in M67 shows that first dredge-up modifies 
surface CN abundances in excess of predictions by 
standard evolution calculations (Gilroy & Brown 
1991; Brown 1987). So almost any amount of mass 
transfer by a giant is likely to modify surface abun- 



dances of blue stragglers formed in this way. Other 
species like Na might also be affected (see Pinson- 
neault 1997, and references therein). The strag- 
gler's companion (the original primary) is likely 
to become either a helium or CO white dwarf as a 
result. 

This is an improvement over lithium and other 
light elements in that it is a pattern change and 
not a destruction of the observable species. If 
the companion's contribution to the spectrum can 
be removed, this could finally provide a means of 
determining the formation path a blue straggler 
took. 

5. Conclusions 

Our spectroscopic analysis of 4 relatively cool 
blue stragglers does not provide definitive evidence 
of composition modifications resulting from the 
formation process. Based on lithium abundance 
patterns as a function of effective temperature for 
main sequence stars in the Hyades and M67, we 
find that current upper limits for the lithium abun- 
dances of all spectroscopically observed blue strag- 
glers are consistent with no mixing during forma- 
tion, even though the upper limits fall below mea- 
sured abundances for turnoff stars. Although the 
observable O lines are formed under NLTE condi- 
tions, the O I triplet equivalent widths for stars in 
our sample and those of Mathys (1991) are consis- 
tent with constant abundance. The three excep- 
tions (S 968, S 975, and S 1434) may be caused 
by emission from or absorption by other sources in 
a binary system. Four of the five blue stragglers 
and all of the main sequence stars have projected 
rotational speeds of less than 20 km s _1 , while the 
fifth straggler (S 975) has vs'mi « 50 km s _1 . 

Although we are unable to determine the ex- 
istence of binary companions to some of the blue 
stragglers in our sample, we are able to provide ar- 
guments about the degree of contamination of the 
spectra of most of the stragglers. We confirm the 
spectroscopic detection of a binary companion to 
the straggler S 1082. From our spectra, we mea- 
sure a projected rotational speed of 90 ±20 km s _1 
for the secondary, and find that its radial velocity 
amplitude is at least 25 km s _1 . The lack of a 
pattern in the radial velocity data (and the vari- 
ability on the order of hours) provides additional 
evidence that this system is currently undergoing 
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mass transfer. The primary is found to be 100 
- 200 K hotter than its companion. S 1082 may 
in the end provide us with definitive proof that 
binary mass transfer can produce blue stragglers. 

We would like to thank the anonymous referee 
for a number of very useful comments that im- 
proved the discussion in this paper. E.L.S. would 
like to thank R. Taam for support (under NSF 
grant AST-9415423) and R. Taam and P. Etzel 
for helpful conversations during the course of this 
work. 
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Fig. 1. — The Fan et al. (1996) data is presented 
here with our program stars marked as open boxes. 
The open triangles are the two stars analyzed by 
Mathys (1991). 
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Fig. 2.— Spectra from JD 2451619 for S 1082 and 
S 984 are shown uncorrected for radial velocity 
variations. A line is drawn at the continuum for 
each spectrum. The three narrow lines are due to 
Fe I. A broad component in S 1082 can be seen 
in each spectrum but not in S 984. The broad 
component of S 1082 moves from blue to red over 
the 5 hours of spectra persented here (bottom to 
top). 
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Fig. 3. — The phased rv curves for the three de- 
tected binaries M67 S 1082, S 984, and S 821. The 
open points represent data from Mathieu et al. 
(1986) while the filled symbols represent veloci- 
ties from this work. Phased curves of three peri- 
ods are presented for S 1082. The horizontal line 
represents the mean cluster velocity (Girard et al. 
1989). 
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Fig. 4. — The combined equivalent widths (EW) 
of the O I 7773 triplet taken from Mathys (1991; 
filled squares) and this study (filled triangles) are 
plotted against B-V color. Lines connect the stars 
in common between the two samples. The large 
crosses arc taken from the Varenne & Monier 
(1999) study of oxygen in Hyades main sequence 
stars. The open and filled circles represent abun- 
dance analyses with constant oxygen abundance 
([O/Fe] = 0.0) for several temperatures using the 
NLTE corrections of Gratton et al. (1999) and 
Faraggiana et al. (1988), respectively. 
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Program Star Radial Velocities 
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Table 3 — Continued 
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Note. — Li EW limits quoted here were taken from Pritchet & Glaspey (1991) 
for 1082 and 997 because their limits were smaller than ours. 
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Table 4 
Program Star Abundances 



star T eff g 
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Li 



[C/Fe] 



[N/Fe] 



[Na/Fe] 



[Mg/Fe] 



[Ca/Fe] 



[Ni/Fe] 



Blue Stragglers 



Our Data: 

Probably Minimal or No Contamination by Companions: 
984 6170 3.9 0.08 ± 0.03 2.20: -0.14 ±0.08 
2204 6650 4.6 -0.05 ± 0.06 <2.00 0.05 ± 0.06 
Unknown Contamination by Companion: 
975 6820 4.4 0.02 ± 0.06 <2.50 -0.09 ± 0.09 
997 6675 4.4 -0.06 ± 0.03 <1.60 -0.03 ± 0.08 
Sure or Probable Contamination by Companion: 
1082 7050 4.5 -0.25 ± 0.05 <1.80 -0.01 ± 0.05 <0.0 
Mathys (1991) Data: 
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Note. — eb — eccentric binary; cb — circularized binary; acc? — accretion occurring in binary?; Am — known Am star 



Table 5 

Program Star Oxygen Abundances 
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-0.04(0.08) 0.16+/-0.10 


-0.10 


0.06 


815 


6275 


4.2 


1.8 


-0.05(0.04) 0.02+/-0.05 
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